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Abstract. Many Lab-on-Chip applications require sample pre-treatment systems. Using electric fields 
to perform cell lysis in bio-MEMS systems has provided a powerful tool which can be integrated into Lab-on-a- 
Chip platforms. The major design considerations for electro-lysis devices include optimal geometry and 
placement of micro-electrodes, cell concentration, flow rates, optimal electric field (e.g. pulsed DC vs. AC), etc. 
To avoid electrolysis of the flowing solution at the exposed electrode surfaces, magnitudes and the applied 
voltages and duration of the DC pulse, or the AC frequency of the AC, have to be optimized for a given 
configuration. Using simulation tools for calculation of electric fields has proved very useful, for exploring 
alternative configurations and operating conditions for achieving electro cell-lysis. To alleviate the problem 
associated with low electric fields within the microfluidics channel and the high voltage demand on the contact 
electrode strips, two “control plates” are added to the microfluidics configuration. The principle of placing the 
two controlling plate-electrodes is based on the electric fields generated by a combined insulator/dielectric 
(glass/water) media. Surface charges are established at the insulator/dielectric interface. This paper discusses the 
effects of this interface charge on the modification of the electric field of the flowing liquid/cell solution. 


1. Introduction 

Many Lab-on-a-Chip applications require sample pre-treatment systems. Using electric fields to 
perform cell lysis in bio-MEMS systems has provided a powerful tool which can be integrated into 

Lab-on-a-Chip platforms [1-3]. The major design considerations for electro-lysis devices include 

optimal geometry and placement of micro-electrodes, cell concentration, flow rates, optimal electric 

field (e.g. pulsed DC vs. AC), etc. The electric field strength must also be moderated to lyse the cell 
but at the same time not destroy cell components that are to be studied. To avoid electrolysis of the 
flowing solution at the exposed electrode surfaces, magnitudes and the applied voltages and duration 
of the DC pulse, or the frequency of the AC, have to be optimized for a given configuration [2]. Using 
simulation tools for calculation of electric fields has proved very useful [2-4] for exploring alternative 
configurations and operating conditions for achieving electro cell-lysis. This paper describes a 
numerical modelling development of a continuous-flow microfluidic device involving using a pair of 
control plates to produce desired electric field for cell lysis. The principle of placing the two 
controlling electrodes is based on the electric fields generated by a combined insulator/dielectric 
(glass/water) media. In fluid flows, both the permittivity and the conductivity affect the flow [5], 



Based on electromagnetic theory [2], surface charges are established at the insulator/dielectric 
interface. The presence of this interface charge will modify the electric field of the flowing liquid/cell 
solution. In the next section, we summarized the leaky dielectric/electric field equations and the 
numerical modelling approach. 

2. Modelling Approach 

2. 1 . Electric Field Equations and Numerical Method 

The governing electric current continuity equation for the lysing device is given as: 

^-fV-J = 0 (1) 

dt 

Where p is the charge density and J is current density. Using the Gauss’ law and introducing the 
electric field strength, E and electric potential, <p , we have 

— (V-eVcp) + V-(aV(p) = 0 (2) 

dt 

Where, s is the permittivity and ct is the conductivity of the media. The electric current density 
is J = oE . 

For numerical solutions, two separate equations are solved in different regions of the calculation 
domain such that: 

V.(-— +e )Vcp n = V • e Vq>°, conductor regions 
• e 0 ( 3 ) 

Vs 0 s r V(p = p, dielectric material regions 


Where, e 0 is the permittivity of vacuum; e r is the relative permittivity with e = s 0 s r ; n indicates the 
solution at the new time step and o indicates the solution at the previous time step. Sub-iterations are 
performed at each time step to calculate the electric potential for the new time step. In the absence of 
an externally specified surface charge density, the conductor-insulator interface condition is a part of 
the solution, enforced using the jump condition: 


3ps 

dt 


int, dielectrics -side 


= V • a' 


int, conductor-side 


(4) 


where p s is the surface charge density at the interface. The implementation for jump condition is the 
total current from equation (4) across the interface: 


J* =(— + £ r )Vcp n -8 r V(p 0 (5) 

£ o 

Then across the interface of conductor/insulator, we have 

J* = (— — ■ + £ )Vcp n -e r V(p° = constant , (6) 

s o 1 

Equation (6) can also be obtained from the first equation in equation (3).This gives us the following 
discretized equation across the interface: 
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Where, subscript c, d, and f mean conductor, insulator and interface. 8n c and 8n d are the normal 


distance from conductor side to interface and insulator side to interface and E° is the normal electric 

field strength at previous time step. It is obvious that the last two terms in equation (7) is the surface 
charge at the interface at the previous time step: 




( 8 ) 


We can solve for interface potential from equation (7): 
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From here, the divergence term of right hand side of equation (3) can be calculated accurately. 

After applying Gaussian theorem at the interface, the charge accumulated at the interface is calculated: 

( 10 ) 

(eE-e 2 E 2 = p s2 

Where, E is electric field strength, p s is surface charge density. 1 refers to the interface with 

conductor 1 and 2 refers to the interface at conductor 2. The above are the basic formulation for the E- 
field calculations in the combined materials. The formulation includes time-domain calculations of the 
electric field and currents in the material. The above numerical method was incorporated into CFD- 
ACE-+- [8] for the electric field calculations. 

2.2 The Microfluidics Cell Lysis Device 

A simple microelectrode geometry as shown in figure 1(a) was the starting point of constructing a 
continuous-flow device. In the figure, also shown in figure 1(b) for the cross-section view, high 
voltage platinum micro-strips are placed 



Figure 1 (a), Basic microfluidic channel with electrodes; (b), Cross section of the electro-lysis channel with thin 

electrodes adjacent to the liquid. 


on the top edge of the channel. Figure 1(b) shows the electric field when 10V DC is applied across the 
two electrodes. It can be seen that at the two lower comers, the field strength reduces to much lower 
level (less than 10 5 V/m). Depending on the cell type, field strength greater than 3x1 0 5 V/m [3, 4, 7] 
should be maintained for cell lysis. 

To alleviate the problem associated with low electric fields within the microfluidics channel and 
the high voltage demand on the contact electrode strips, two “control plates” are added to the 
microfluidics configuration. These two plates are placed at about 500 micro meters above and below 
the channel. The theoretical background for placing the two controlling electrodes is based on the 
electric fields generated by a combined insulator/dielectric (glass/water) media. In fluids, both the 
permittivity and the conductivity affect the flow [5], Based on electromagnetic theory [6], surface 
charges are established, at the insulator/dielectric interface. Based on the numerical modelling results, 
this surface charge and the electrical double layer established will generate sufficient electric field 
(perpendicular to the electric field generated by the lysing electrodes) at the comer regions. The 
suggested operating conditions summarized below are based on the conductivity of the dielectric 
media (water), as well as relative permittivity of both insulator (glass/Si02) and dielectric material 
(water). Several operating methods are possible. For example, either AC or square pulses [1,2] can be 
applied at the electrodes. To avoid water electrolysis [2], AC was used in this analysis. Due to the 
large difference in time constants (estimated from the ration of permittivity and conductivity), a DC 
was applied across the control plates. Here we present the results of AC of 20V across the electrodes 
and a constant 400V across the control plates. The AC was simulated as a cosine wave with 0.1MHz. 
In figure 2, the electric field magnitudes are shown at the beginning of the cycle (2.a) and at the half 
cycle (2.b). It can be observed that the low-strength areas at the two lower comers have been much 
reduced. In figure 3, the surface charges at both top and bottom channel walls are plotted at the 
beginning, at 14 , and at half of the AC cycle. 


3. Concluding Remarks 


In this study, we have utilized a numerical simulation tool to investigate the design of a continuous 
microfluidics based cell lysis device. Original lysing channel configuration with two thin electrode 
strips wall-flushed at both upper portion of the channel creates low electric field outside of “kill zone” 
[4] at lower portion of the channel. The proposed new configuration by incorporating a pair of control 
plates placed above and below the channel has shown its effectiveness in increasing the electric field 
strength. Comparison with experimental testing is currently underway to verify the proposed operating 
conditions. 




Figure 2. Cell lysis channel with control panels (two panels above and below the channel not shown); (a) electric 
field magnetite at the beginning of the cycle, (b) at half cycle. 


-1.SE44 
-LIE 44 

N 

4 

£ -2.0E -04 

<5 

-2.X 44 
-2.6E44 

-2JE4f LDE-* 4.0E44 &QE-K l.*E4# 2.*E4t 



X,m 



-LIE44 
-22E44 
^-2j6E44 
-3JE44 

-15E-05 -L5E05 -54EI6 5jQE06 L5E45 2.5E45 

-LIE-14 
-UE-04 

CM 

S -2.2E-04 
|.24M4 
-3.4E-I4 

-4JE-04 

2.50E45 -L50E-O5 5J0OE46 5.00E06 1.50E85 2.50E 

X, » 






Figure 3. Surface Charges at the top wall (left column) and bottom wall (right column) of the microchannel; (a) 
at the beginning of the cycle (top row), (b) at quarter cycle (middle row), (c) at half cycle. 
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